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Abstract. Over the past decade we have seen an explosion of demonstrations 
of quantum coherence in atomic, optical, and condensed matter systems. These 
developments have placed a new emphasis on the production of robust and optimal 
quantum control techniques in the presence of environmental noise. We discuss the use 
of dynamical decoupling as a form of open-loop quantum control capable of suppressing 
the effects of dephasing in quantum coherent systems. We introduce the concept of 
dynamical decoupling pulse-sequence construction as a filter-design problem, making 
connections with filter design from control theory and electrical engineering in the 
analysis of pulse-sequence performance for the preservation of the phase degree of 
freedom in a quantum superposition. A detailed mathematical description of how 
dephasing and the suppression of dephasing can be reduced to a linear control problem 
is provided, and used as motivation and context for studies of the filtration properties of 
various dynamical decoupling sequences. Our work then takes this practical perspective 
in addressing both "standard" sequences derived from nuclear magnetic resonance 
and novel optimized sequences developed in the context of quantum information. 
Additionally, we review new techniques for the numerical construction of optimized 
pulse sequences using the filter-design perspective. We show how the filter-design 
perspective permits concise comparisons of the relative capabilities of these sequences 
and reveals the physics underlying their functionality. The use of this new analytical 
framework allows us to derive new insights into the performance of these sequences 
and reveals important limiting issues, such as the effect of digital clocking on optimized 
sequence performance. 
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1. Introduction 

The needs of the field of quantum science and quantum information have evolved 
over the past decade. Just ten years ago the bulk of theoretical and experimental 
effort in the community was focused on the realization of new quantum coherent 
systems using condensed-matter devices, optics, and atoms. The successes of the 
field [U[2l[3l[ll[5l[6l[Tl[8l[9l HOj HH [121 HH1 HH US] in this short timeframe 
have led to the emergence of a growing need for error-resilient quantum control 
techniques. Such techniques are enabling for probing studies of the dynamical evolution 
of quantum systems, the realization of high-fidelity quantum logic operations in quantum 
information, and the production of useful quantum coherent technologies. 

The "continuous" information representation used in quantum technologies means 
that small deviations in a quantum system's state from the nominally ideal state yields 
an increased probability of error in a particular basis [18]. This makes quantum states 
particularly sensitive to unwanted and unavoidable environmental couplings, resulting 
in state randomization through the process of decoherence. Suppressing decoherence 
has thus become a primary focus of the community as interest shifts away from basic 
demonstrations of quantum coherence towards attempts at realizing useful quantum 
coherent technologies [191 I2Q] . 

The field of quantum control [21] addresses a broad set of challenges pertaining 
to the robust and efficient manipulation of quantum systems and processes. So- 
called "open-loop" control methods, in particular, rely on feedback-free Hamiltonian 
engineering through the application of time-dependent control sequences, and leverage 
strengths from longstanding fields such as optimal control theory and nuclear magnetic 
resonance (NMR) [22J. The latter holds particular importance as modern magnetic 
resonance imaging (MRI), a biomedical application of NMR, is largely based on pulsed- 
control techniques designed to overcome the semiclassical analogue of decoherence in 
a macroscopic ensemble of nuclear spins (for instance in a tissue sample). The many 
lessons learned in this field, beginning with Hahn's discovery of the spin echo in 1950 [23J, 
can be brought to bear in the challenge of realizing error-resistant hardware for quantum 
computing applications. 

Beginning in the late 1990s significant interest developed in the use of modified 
concepts from NMR and open-loop control theory for the suppression of decoherence 
and error in quantum computing applications [2H [2S1 [23 EH1 E01 EH E21 E31 [M] ■ 
These protocols, known as dynamical decoupling (DD), are extensions of Hahn's spin 
echo [221 ES] , modified to suit general decoherence processes in qubit systems including 
the influences of a bath that itself exhibits quantum dynamics. In these techniques a 
prescribed sequence of control operations at the physical level enables environmental 
fluctuations to be coherently averaged out, thus preventing errors from accumulating 
and preserving quantum information. The primary degrees of freedom in this approach 
to quantum error suppression are the number of control operations, n, and their relative 
timing, leading to a wide range of possible DD sequence constructions. 
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Altogether, the true power of dynamical error suppression techniques is derived 
from the facts that: 1) only unitary control operations are employed - no measurement 
capabilities nor encoding overheads are required; 2) errors may be suppressed regardless 
of magnitude, as long as appropriate time-scale separations are obeyed; 3) sequences 
may be optimized for a wide class of noise processes. These facts imply that a sequence 
of control operations may be applied without an experimentalist having a detailed, 
real-time quantitative knowledge of the ambient decohering environment, and if applied 
correctly, an arbitrary quantum state of interest will be preserved. 

Using DD it is thus possible to delay the decay of quantum information and suppress 
decoherence-induced errors to the levels needed for large-scale quantum computation or 
other applications of quantum systems. The qubit "refresh" techniques to be studied 
in this work are similar to firmware refreshes of DRAM cells used to compensate for 
classical information leakage [37] . If quantum information is protected from decoherence 
in a similar way using dynamical error suppression protocols in support of some ultimate 
computation or other application, it is appropriate to consider these protocols a form 
of universally applicable "Quantum Firmware." 

In this manuscript we focus on a presentation of the analysis of dynamical 
decoupling pulse sequences as a problem in filter design, in which we map familiar 
concepts from electrical engineering such as the 3 dB point, high-pass/band-pass 
filtering, and low-frequency rolloff (for a high-pass filter) to performance evaluation 
of various sequences for the problem of pure dephasing. The objective is to give 
an accounting of state-of-the-art techniques in dynamical decoupling, but rephrasing 
the problem in a new framework which has an accessible and practical interpretation, 
including a quantitative presentation in the context of control theory. Our work builds 
on several recent theoretical efforts aiming to explain decoherence and decoherence 
suppression in terms of spectral functions [2S1 ESI EE]- Further, it complements 
alternative approaches that have been presented recently, using the analogy of optical 
interference from a diffraction grating to describe dynamical decoupling sequence 
performance [3D] . 

Using the filter-design perspective we study familiar pulse sequences such as the 
Carr-Purcell-Meiboom-Gill (CPMG) multipulse spin echo, and Periodic Dynamical 
Decoupling (PDD) pU [25] and compare against recently discovered sequences such as 
Uhrig Dynamical Decoupling (UDD) [38], Locally Optimized Dynamical Decoupling 
(LODD) pE], Optimized Noise Filtration by Dynamical Decoupling (OFDD) [ID] . 
Bandwidth- Adapted Dynamical Decoupling (BADD) [H], and others. In all cases 
we leverage insights from filter design to concisely explain the performance of the 
various sequences when applied to the mitigation of decoherence. Our perspective also 
permits us to extract new information on the effects of timing constraints on sequence 
performance, exclusively through study of the modification of the sequence filters. 

Dynamical decoupling can be used to preserve arbitrary quantum states in the 
presence of general decoherence processes - leading to depolarization. In this work we 
focus exclusively on the case of pure dephasing, which while constrained, is of great 
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importance to the community. Further, we do not address concatenated sequences [30] 
or their performance characteristics in suppressing multi-axis decoherence, and instead 
focus on the construction of pulse sequences that could form a single concatenation level 
in a more complex sequence [121 H3j. The applicability of the filter-design concepts to 
more general forms of decoherence remains a subject for future study but is supported 
as a general proposition by past theoretical work [29|. 

Quantum systems come in many forms, but we will focus on the quantum bit, or 
qubit, that is most prevalent in quantum information science [18J. This abstraction 
of a quantum mechanical two-level system is useful for information processing, but 
also serves as a model system for understanding the dynamical evolution of quantum 
coherent states exposed to noisy environments. We will frequently employ discussions 
of technical details derived from quantum information processing as it represents one of 
the best-known applications for quantum technologies to date. 

The remainder of this manuscript is structured as follows. In Section [2] we review 
a mathematical formalism for qubit dephasing that is generalized to include the effects 
of a dynamical decoupling pulse sequence. This is followed by Sec. [3] which forms 
the core theoretical section of the manuscript. Here, we describe how dephasing can 
be reduced to a problem in linear control, using insights from control theory, and 
explain how dynamical decoupling pulse sequences behave as filters in the task of 
noise suppression. The frequency-domain performance of a given dynamical decoupling 
sequence is presented in |4j and we compare the performance of various pulse sequences 
and discuss relevant operability regimes of interest. Using these insights, in Sec. [5] 
we present new analyses of the effects of timing constraints on filter characteristics, 
including a detailed study of the incorporation of physically realistic (noninstantaneous) 
control pulses and discretization of time, for instance due to digital clocking. Sec. [6] 
reviews techniques to create pulse sequences with spectral filtering characteristics 
optimized for a given noise environment using search algorithms. The manuscript 
concludes with a summary and brief future outlook. 

2. Error Accumulation Due to Environmental Noise 

One of the most pernicious error sources in quantum coherent systems is the 
phenomenon of decoherence, in which random qubit errors accumulate due to unwanted 
environmental coupling. In this section we will describe the effects of decoherence on 
a qubit system quantitatively, in order to provide physical insight to the dynamical 
evolution of a quantum state. 

The effect of decoherence can be seen by examining the temporal evolution of a 
qubit initially in a superposition state, e.g. 



where \i) are the qubit basis states and are properly normalized, complex 
probability amplitudes that permit a more compact representation of the state. In this 




(1) 
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picture an arbitrary qubit state in this two-dimensional Hilbert space can be represented 
as a vector on a sphere of unit radius. This so-called Bloch-sphere representation (lower 
panel, Fig. [I]) illustrates how both the relative probability amplitudes and the relative 
phase between the qubit basis states uniquely define a particular superposition state. 
Unwanted environmental coupling has the effect of randomizing the two variables in 
which state information is encoded, Cf and c±. 

One may consider two general classes of decoherence processes: Longitudinal energy 
relaxation and transverse dephasing. In the former process, the probability amplitudes 
of the qubit state are affected (i.e. 6 or |cj| are randomized), forcing the qubit state to 
move along a meridian of the Bloch sphere. The characteristic time over which a two- 
level system undergoes energy relaxation is known as T\ [35J, borrowing from the NMR 
literature. Transverse dephasing involves randomization of the relative phase between 
the basis states. In an ensemble average - for instance, averaging over many different 
experimental outcomes - such phase randomization leads to a decay of coherence in a 
characteristic time T^. In total, these two decoherence processes limit the useful lifetime 
of the system to T 2 _1 = T^ 1 + {2T l y l . 

The dynamics of a qubit in a general decohering environment may be represented 
in a semiclassical picture are governed by a Hamiltonian of the form 



where cr* represent the Pauli matrices, HQ is the unperturbed qubit energy splitting, 
and a, (3 represent random fields imparted by the environment. It is generally sufficient 
to consider only the transverse and one longitudinal process due to the noncommuting 
nature of the Pauli operators, hence the omission of a term dependent upon ay. 

Application of the time-evolution operator for time t with a = (5 = leads to a 
phase accumulation Qt/2 in the Schrodinger picture. This phase evolution must always 
be tracked in a quantum control setting; any arbitrary qubit state accumulates phase 
due to this term in the Hamiltonian. It is often convenient, then to write 



where the initial phase of the qubit state at time t = is absorbed into the Cf^. In 
this way the phase evolution due to Schrodinger's equation is captured explicitly in the 
variable </>. 

From the perspective of an experimentalist interested in quantum control, qubit 
coherence is often dominated by transverse dephasing processes. In a pure dephasing 
picture a = 0, f3(t) ^ 0, and the state 1*1/(0) evolves to 



where the phase dependent upon f2 has been dropped by transforming to the rotating 
frame. This is represented in the density matrix formalism 



H = a{t)a x + [- + f3{t)]cr z , 



(2) 



¥(*)) = ct |t) e^ 2 + c^ 2 ||) 



(3) 



(4) 



Pit) 




2i J* P(t)dt 




(5) 
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where the off-diagonal matrix elements give the system coherences, which decay to zero 
in an ensemble average: decoherence. For the remainder of this manuscript we shall 
consider exclusively the case of pure dephasing. 

Classical noise in common control or environmental parameters forms a dominant 
source of dephasing in many quantum control settings; typical dephasing mechanisms 
include magnetic field fluctuations for atomic systems [I] , charge fluctuations for charge- 
qubits in solids [44J , and the effective Overhauser field arising from randomly fluctuating 
nuclear spins in semiconductor spin qubits [45] 06]. Experimentalists must also consider, 
however, the degree to which one may experimentally realize the implicit transformation 
to the rotating frame. In a physical sense this transformation corresponds to always 
keeping track of the phase of the qubit relative to a master oscillator. We see, therefore, 
that aside from fluctuations in external parameters, noise in the master oscillator [THUS] 
can also contribute to an effective dephasing between the control hardware and the 
qubits under test, and that this dephasing may be represented as a contribution of 

pa) m. 

Using these insights we focus on a semiclassical representation of the influence of 
noise, rather than a full quantum mechanical treatment as may be found in Ref. [49, 34J. 
The ensemble-averaged phase accumulated between the qubit basis states due to random 
fluctuations in /3(t) may therefore be equivalently represented as a time-integral over 
(3{t) or as a frequency integral over the power spectral density. 

We consider a formulation for measuring the coherence for a dephasing Hamiltonian 
presented originally by Uhrig [381 EQ] and Cywinski [51]. Starting with a state |^ ) 
oriented along the F-axis in the equatorial plane of the Bloch sphere, the state 
accumulates a random phase due to environmental noise. The coherence of this state 
after time r is 

W(t) = |K)(t)I = e~* (T) (6) 

where angled brackets indicate a quantum-mechanical expectation value and the overline 
indicates an ensemble average. 

Analysis may be conducted in either the time-domain or the frequency domain. It is 
often convenient to transform the time-domain noise term to the frequency-domain 
power spectral density, 

/oo 
e*" (/3(t + T)P(t))dT, (7) 
-oo 

in which case we may express 

X(r) = 2 - J^F(ujr)du. (8) 



Here F(ut) describes the spectral dependences of the experiment being performed and 
is known as the "Filter Function." 
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2.1. Free- Induction Decay 

The first kind of experiment of interest involves allowing a state to freely precess under 
the influence of a dephasing Hamiltonian, allowing one to effectively probe the temporal 
dynamics of the state as in Eq. |4} After some time r, the state accumulates a random 
phase due to the time-integral of /3(t). In an ensemble average this produces a decay of 
coherence, similar to the free-induction decay (FID) experiments of NMR. Here, we are 
interested in the total average accumulated phase-error, obtained by integrating over 
the experimental duration or over the relevant window in frequency space, and thus are 
only marginally concerned with the spectral content of Sp(oj). It has been shown that 
F^ fid \ujt) = sin 2 (cot/ '2) j36j [50], and the coherence integral may then be written 

X(r)= 2 /^sin>r/2)du,. (9) 

7T J CO 


In the limit of low-frequency-dominated fluctuations this expression gives a familiar 
Gaussian decay envelope to the coherence [53]. A primary motivation in the following 
discussion is to extend this coherent lifetime and suppress decoherence-induced errors. 



2.2. Coherence with Dynamical Decoupling 

Dynamical decoupling involves the sporadic application of control pulses in order to 
average away the effects of coupling between the quantum system of interest and the 
environment. We focus on an approach that utilizes the quasi-periodic application of 
parity-reversing pulses to effectively time-reverse the accumulation of random phase 
in successive free-precession periods. More general approaches, including randomized 
decoupling [54], continuous arbitrary modulation [29j 155]. or other Pauli rotations [56] 
may also be considered, but analyzing their construction in the context of filter design 
will be the subject of future work. 

Assume arbitrary dephasing-noise as described above and a sequence of n parity- 
reversing 7r-pulses (Ily = e*5' Jy = iay) during which the instantaneous pulses are applied 
at fractions 5j, j = 1,2, ...n, of the total free-evolution time r. In this case the time 
evolution of an initial state is given by 

|*(t)> = e-^^^^^^.-nye-^C^^^nye-^C^ 4 ^*' |*( )> .(10) 

In the absence of 7r-pulses the total random phase that the qubit accumulates is given 
by 4>{t) = Jq [3(t)dt as described in the previous section. After commuting every other 



operator Ily to the left once in Eq. (10), we see that due to the 7r-pulses this random 
phase becomes instead 

0(r) = (-l) n / f3(t')dt'...+ (3(t')dt'- f3(t')dt'. (11) 

In this construction, the state freely evolves, accumulating a random phase, a 7r-pulse is 
applied and the system continues to evolve again accumulating phase, but with opposite 
sign to that in the first period, and so on. For an arbitrary time sequence of j3(t), 
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Figure 1. Schematic representation of a single-axis dynamical decoupling pulse 
sequence and its relation to the filter function, (a) Schematic representation of the 
time-domain application of tt pulses. Pulse width is represented as being nonzero for 
clarity, and are applied about the Y-axis. Total sequence duration is r. (b) Graphical 
representation of the toggling frame for phase accumulation corresponding to the 



integrals in Eq. 11 Phase accumulation "winds" and "unwinds" during successive free- 
precession periods. Frequency information about the performance of a given dynamical 
decoupling pulse sequence is captured in the temporal form of y n (t). 



t G [0, r] the phase evolution of |^) due to the time integral of (3(t) is modulated by a 
"sampling function," y n (t), which takes alternating values ±1, with transitions occurring 
at the temporal-location of the 7r-pulses. Dynamical error suppression therefore involves 
breaking up the phase-evolution into a series of shorter free-evolution periods in which 
phase sequentially and repeatedly "winds up" and "unwinds" (Fig. [Tp, producing an 
effective averaging of the phase accumulation through the sampling function. 

In the case that /3(t) is static but nonzero over the relevant interval [0, r], a single 
7r-pulse located at time r/2 can perfectly cancel phase accumulation, constituting the 
"spin-echo" as developed originally by Hahn. Any variation of (3(t) between the two free- 
precession periods reduces the net coherence of the system as the phase accumulation 
will be slightly different in the two periods. It is known that by applying a series 
of n 7r-pulses in a particular, evenly spaced sequence, one may form the Carr-Purcell 
multipulse spin echo (or CPMG given a particular phase relationship between the applied 
it pulses and the qubit state). Here the phase accumulation is broken up such that in 
successive free-evolution periods f3(t) is kept nearly constant, thus maintaining phase- 
coherence. Again, variation of /3(t) between successive free-evolution periods reduces 
the net coherence of the qubit at the conclusion of the sequence. 

Unlike an FID experiment, the performance of a dynamical decoupling pulse 
sequence is extremely sensitive to the frequency spectrum of 0(t). An applied dynamical 
decoupling sequence temporally modulates a qubit's phase accumulation, as presented 
above; the total free-precession time period is split into smaller time bins separated by 
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7r-pulses. Slow fluctuations in f3{t) are corrected via the applied 7r-pulses, while more 
rapid fluctuations between free-precession periods contribute to the net dephasing of the 
qubit. The relevant frequency components are set by the total sequence length r and 
the interpulse periods. As we will see, it will be extremely useful to consider dynamical 
decoupling pulse sequences as spectral filters. 

3. Sequence Construction as Filter Design 

The multipulse spin echo (to which we interchangeably refer as CPMG) is a simple 
dynamical decoupling pulse sequence with a long history of use in NMR. However it can 
face significant performance limitations arising from its construction. As discussed in 
the last section, the presence of noise that fluctuates rapidly relative to the interpulse 
spacing diminishes the efficacy of the multipulse spin echo. In quantum information 
experiments we aim to maximally preserve coherence rather than simply extend the 
lifetime of a signal, as in NMR, in the presence of noise characterized by arbitrary 
power spectral density. Our present task will be to use the degrees of freedom available 
to us in order to create new pulse sequences with augmented performance for the specific 
case of noise which is not well approximated as quasi-static. We will accomplish this by 
considering pulse sequence construction as a problem in filter design. 

3.1. Filter Basics 

Electrical engineering and digital signal analysis tasks often require an input signal to 
be mapped to an output signal through a known transformation. This is frequently 
enacted in one variable and analyzed in its conjugate. A prototypical example is the 
action of changing the spectral characteristics of a time-domain signal, for instance, by 
removing high-frequency components above a given cutoff. This approach is generally 
useful in the removal of extrinsic interference in a signal processing task - for instance, 
removing mains power frequency modulation of an electrical signal by using a low-pass 
filter. In analog electronics, filtering may be achieved using hardware components - for 
instance networks of passive devices such as resistors and capacitors. However, given a 
mathematical representation of a time-signal, a Fourier-domain mathematical filter may 
also be applied in the case of digital signal processing. For instance, one may perform 
an FFT on a time-domain signal, remove certain spectral components in software (via 
the action of the filter), and transform back to the time-domain. 

It is instructive to introduce a few basic terms from filter design [57] in order to 
make plain the corresponding features of classical filtering and the action of dynamical 
decoupling sequences. Here, we present a brief description of salient characteristics that 
is by no means exhaustive. 

• The Transfer Function: this describes how a given input signal is transformed to 
an output signal by the action of the filter. For a generalized, arbitrary filter we ex- 
press the application of the filter as the convolution of the input signal with the time 
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domain transfer function of the filter. This is more commonly expressed through 
the Laplace transform as y ou t{ s ) — G(s)y in (s), where s is a complex variable, and 
G(s) is the transfer function in the fourier domain. In discrete time this may be 
represented in matrix form. 

• Frequency Response/ Filter Gain: used to characterize the aggregate behavior of 
the Transfer function to a given impulse as expressed in the frequency domain, and 
describes how particular frequency components of an impulse are modified with ei- 
ther positive (enhancement) or negative (suppression) gain. For instance, this may 
include high-pass, low-pass, band-pass, and notch behavior. 

• Pass-band and Stop-band: the regions of frequency space corresponding to trans- 
mission and attenuation of relevant frequency components by the filter. 

• Filter Order: a mathematical parameter used in the construction of the filter which 
often corresponds to the number of repeated physical units in, e.g. an electrical 
analog filter. The filter gain is often conveniently expressed in terms of the filter 
order. 

• Roll- off: describes the steepness of suppression of a given frequency range below 
(above) cutoff in a high-pass (low-pass) filter. Steeper roll-off corresponds to more 
efficient suppression of noise in the stop-band. An 77-order filter has 61] dB/octave 
= 2O77 dB/decade suppression of signals below (above) cutoff. 

• 3-dB Point: frequency below (above) cutoff at which frequencies are suppressed 
by 3 dB relative to the pass-band in a high-pass (low-pass) filter. For practical 
purposes this is often associated with the cutoff frequency. 

• Ripple: nonuniformities in the filter gain lead to frequencies being transmitted in 
the pass-band with varying amplitudes. 

To understand how the concept of filter design may be applied to dynamical 
decoupling and quantum control, we must first establish a correspondence between 
the mathematics of filtering (from the perspective of control theory) and dynamical 
error suppression. This requires an explanation of how dephasing due to random noise, 
characterized by a power spectral density may be reduced to a linear control problem. 
The following subsection will use terminology derived from the control theory literature 
in order to demonstrate the appropriate relationships. We will then return to the 
language used more commonly in the physics literature. 
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3.2. Dephasing and DD as linear control problems 

The linearity of the dephasing of the qubit may be understood as follows. We saw in 
the previous section that at the end of a dynamical decoupling sequence that starts at 
t — t and finishes at t we have, before averaging over the noise, 

(oy) = Re{exp(-2i0(t))} = Re jexp [^ 2i f ( ) M* ~ t')/3(t')dt'^j J . (12) 

Here h is the impulse response function implicitly defined in equation (11), which defines 
the filter that is applied to the noise to produce the phase 0, and is equivalent to 
the sampling function y n (t) for a chosen pulse sequence. Writing things in this way 
emphasizes that there is a linear dynamical map from the noise /3 to the phase of the 
coherent qubit oscillation, and suggests the use of techniques from filtering and control. 
It is important to note though that it is not possible to generate a general impulse 
response h from the sorts of dynamical decoupling pulses given in equation (11) so the 
class of possible linear evolutions is constrained. 

When we average over the noise we find our figure of merit 

W = e-* (r) = e -^ 2(T)) (13) 

and it is clear that we can optimize our pulse sequence of length r if we minimize 
(0 2 (r)). We regard this as our objective function to be minimized. Intuitively h should 
be designed so that (f> does not respond sensitively to (3. In control theory language, 
is not very controllable by the noise input /3. 

In control theory there is a standard description of this sort of problem as 
follows [58J : consider a physical system with states described by a vector x undergoing 
linear dynamics and driven by some control or noise signals given by u that produces 
some output y that is a linear function of x 

x = Ax + Bu{t) (14) 

y = Cx, (15) 

which may be solved to give 

x(t) = e At f e- At ' Bu(t')dt' (16) 
J o 

y (t) = Ce At f e- At ' Bu(t')dt' . (17) 
Jo 

Applying the Laplace transform we find y(s) = G(s)u(s), where G(s) = C (si — A)^ 1 B 
is the transfer function for the control operation of interest, / is the identity, and s is a 
complex number. 

In our system we simply have one input u(t) = (3{t) and one output y(£) = x(£) = 
4>(t). We describe our dynamical decoupling experiment as follows. We consider our 
state to be initialized at time t = —oo and evolve unperturbed until time t = — r. At 
this time the Hamiltonian term /3(t) = Bu(t) is turned on and the state allowed to 
evolve until time t = 0, at which point the state evolution is said to be complete. In a 
free evolution process we set G(s) = /, and the state evolves under the application of 
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the dephasing Hamiltonian. 

When considering the prevention of dephasing we are dealing with a linear control 
problem in which we are attempting to implement noise suppression. Our system 
experiences a disturbance defined by /3(t) = Bu(t), and our task, in applying a 
dynamical decoupling pulse sequence, is to filter out the effect of the noise, i.e. to 
reduce the system's controllability by the environment. As a result of our objective 
function we do not care about values of the output for t < 0, only the phase at the end 
of the pulse sequence matters. As a result we are interested in the map from the input 
noise signal u(t) for t < to the state at t — 0, x(0). This mapping is linear, and in 
control theory is known as the controllability operator. We write 

x(0)=<*f[u(t)], (18) 

where we can think of the controllability operator ^ [u(t)] as an asymmetric matrix 
that represents the transform of u under the action of transfer function G(s). The 
expectation value of interest for our control problem is (x T (0)x(0)^ = (u T( t^ T( t^uJ, 
which we seek to minimize. (The angle brackets here indicate averaging over realisations 
of the noise.) This formula is analagous to the construction of time evolution under a 
dephasing environment and an open-loop control protocol to be represented by the 
quantum mechanical operator 0, where the expectation value is represented at time 
t = as (\P(0)| |\P(0)) = Tr (0 |^(0)) (^(0)|). This suggests some further simplifcation 
of the formula for the objective function. We may write, 



(x T x) = (u T ^u) (19) 

= (Tr (? T ?uu T )) (20) 

= Tr(V T ^(uu T )) (21) 

= Tr (V V^ r ) . (22) 



Here the quantity V denotes the covariance matrix of u, playing a role directly 
comparable to Sp(u). In the case of Gaussian white noise this reduces to the identity, 
however this quantity is nontrivial in the case of non-Markovian processes which exhibit 
colored noise (e.g. 1/w). The quantity ^V^ 5 " is in general a matrix, known to 
control theorists as the weighted controllability Grammian, and provides all information 
about how noise with a colored spectral representation, under a control procedure (here 
dynamical decoupling) maps x(— oo) — > x(0). (Weighted because of the effect of V 
to account for the spectrum of the noise). For our class of problems we would wish 
to minimize the eigenvalues of the controllability Grammian by choosing the transfer 
function G such that the system does not respond strongly to the noise. In our case 
the weighted controllability Grammian is simply the number x- 
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3.3. DD in the Frequency Domain 

The control-theoretic task in a dynamical decoupling experiment using continuous 
(discrete) time can be reduced to construction of a F(ujt) (G, or < &) that maximally 
suppresses noise characterized by Sp(u>) (V). 

The time-domain modulation due to a decoupling sequence enacts the desired 
noise suppression through the convolution of (5{t) and y n ,(r). Hence we look to 
the modulation's Fourier transform to provide the relevant spectral information. In 
particular we may write the filter function [381 EH] , 

n 

F(ujt) = \y n (tur)\ 2 = |1 + (-1)™+V" T + 2 £(-l)V^ T | 2 (23) 

This quantity provides all information about how an arbitrary pulse sequence will 
suppress phase accumulation as a function of frequency, and exists in analogy to the 
filter gain. The filter function enters our expression for coherence, W(r) = e~ x< ^ T \ as 
described in Eq. [8] Since the filter function appears as a multiplicative factor of Sp(uj), 
small values of F(ut) in the dominant spectral regime will lead to small values of x( r )) 
and hence coherence W(t) m 1. We see from this presentation that the control-theoretic 
< ^ 72 oc \F{ut) from the physics literature. 

Examination of Eq. [23] indicates that changing n or modifying the fractional pulse 
locations 5j, in an n-pulse sequence can alter its frequency response. The fact that Sj 
is an easily controlled parameter has produced a recent focus in the community on the 
construction of new DD sequences tailored to suppress dephasing noise with arbitrary 
spectral characteristics by modifying the form of the filter gain. 

A sequence design methodology based on modification of the filter function has 
provided a new avenue for the construction of experimentally relevant sequences and 
has allowed a simple physical insight into the action of dynamical decoupling broadly. 
An experimental operator may construct a dynamical decoupling sequence to provide a 
F(ut) that is most useful for the noise present in the experiment - one that minimizes 
the objective function defined above. 

This approach is taken routinely when an experimentalist attempts to improve 
qubit coherence by adding pulses to a CPMG sequence; doing so reduces the minimum 
pulse separation for a given total experimental duration and hence allows correction 
for higher- frequency (more rapidly changing) components of Sp(u). Recent theoretical 
investigations [301 EH EH [591 EI] have demonstrated that substantial improvements in 
achievable memory fidelities in the presence of rapidly fluctuating system baths (and 
sharp high-frequency cutoffs) are possible by resorting to more sophisticated sequence 
design, as will be described below. 

Equivalent to the control theory perspective above, the frequency- domain filter 
function shows how an applied dynamical decoupling pulse sequence behaves explicitly 
as a spectral filter, with gain characteristics captured in the mathematical form of 
F(ut). For low frequencies the phase-factors in the filter function construction are 
small compared to 7r/2 and the terms add destructively to produce the stop-band. 
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In this regime the filter function remains small and noise contributions to x{t) are 
minimized. For the highest frequencies the phase-factors can differ by 7i~, permitting 
constructive interference and yielding rapid fluctuations about the mean value An + 2. 
Frequency components of Sp(cu) in this regime may contribute substantially to x(t), 
yielding dephasing and primarily impacting T 2 . In this way we may think of a dynamical 
decoupling sequence as behaving like a high-pass filter. The filter function in the 
pass-band exhibits significant ripple and can take values greater than unity, producing 
effective enhancement of dephasing due to noise at certain frequencies. 

We define the analog of a classical electrical filter's 3-dB point as the crossover 
between the two regimes described above; it is the value of cut (or u for a given r) giving 
a value of the filter function F n ~ 1, denoted ojpi- This corresponds approximately to 
the "dynamical decoupling limit," |4"Tj 



where r m j n is the minimum interpulse period. This says that dynamical decoupling will 
successfully permit rephasing so long as the highest-frequency component of Sp(u) is 
slow relative to the smallest interpulse period. Returning again to the filter function 
formalism, if Sp(u) has significant components u > upi, the system will be largely 
dephased as the applied sequence cannot decouple the qubit from the rapidly fluctuating 
environment, i.e. the filter gain passes these frequencies approximately unimpeded. 

In many experimental studies of coherence, r is varied and the system's coherence is 
measured. Because r enters the expression for the filter function and sets the frequency- 
response characteristics of a given dynamical decoupling sequence (along with 5j and n), 
variation in r will change the frequency response of the sequence. Thus the accumulation 
of error will not necessarily result in a simple exponential or even Gaussian decay 
envelope in r when dynamical decoupling pulses are applied. Instead, the form of the 
decay can be highly irregular for a particular sequence even showing "revivals" where the 
measured error occasionally decreases with increasing r for certain ranges, depending 
upon the characteristics of the noise [3]. 

4. DD Filter Characteristics 

4-1- Sequences of Interest 

Early DD schemes have largely relied on the simple periodic repetition of 
(approximately) instantaneous pulses. The best known among these are the CPMG 
(after Carr Purcell, Meiboom, and Gill, and referenced above) and Periodic Dynamical 
Decoupling (PDD) sequences. 

The CPMG sequence was originally developed in the context of NMR systems where 
inhomogeneous broadening required an ability to refocus the ensemble Bloch vector as 
it spread out during free precession [22]. However, the same sequence is quite effective 
at suppressing phase randomization due to homogeneous effects (e.g. (3(t)) when noise 




1 



(24) 
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processes are dominated by low-frequency components (e.g. Sp(cu) oc 1/w). An n-pulse 
CPMG sequence is constructed with fractional pulse locations 5j = (j — 1/2) /n (Fig. [l^,) 
wherein the first and last free-precession periods are half the duration of the interpulse 
periods, producing effective "refocusing" of the Bloch vector at the conclusion of the 
sequence. 

Periodic Dynamical Decoupling involves the repetitive application of uniformly 
spaced tt pulses with fractional locations Sj = j/ (n + 1) for n pulses [24J. The sequence 
does not provide efficient refocusing of the Bloch vector, but does perform averaging of 
environment-induced phase accumulation due to its sampling function [51]. However, 
this sequence has the benefit of providing suppression of general decoherence, if properly 
constructed. 

Uhrig first showed that manipulation of the relative pulse locations, Sj, for fixed n 
and r leads to modification of F(ur), providing the ability to tailor the filter function. In 
particular, he analytically derived an n-pulse sequence in which the first n derivatives of 
y n (ur) vanish for wr = 0. Expressed alternatively, given an expansion of j3(t) in powers 
of t, (3(t) = (3o + (5\t + /3 2 t 2 + f3st 3 + the UDD sequence is constructed to suppress the 
first n orders of the expansion, given n pulses, while the CPMG sequence cancels exactly 
only the zeroth and first-order terms for all n > 2 [60, 9j. The resulting sequence, UDD, 
has n pulse locations determined analytically as Sj = sin 2 [7rj/(2n + 2)]. 

The construction presented above allowed Uhrig to tailor the filter function such 
that it provided strong suppression of phase accumulation when noise environments 
possessed significant high-frequency contributions — a dramatic advance over CPMG. 
Uhrig specifically showed [50] that in noise spectra including high-frequency components 
and a sharp high-frequency cutoff, ujj, such as the Ohmic spectrum (Sp(u) oc ojQ(ujd — 
(J)) that may be consistent with a spin-boson model [HI EH] [29], the UDD sequence 
would yield significant gains in performance relative to CPMG [3B1 EO]- By contrast, 
in the presence of noise with a soft high-frequency cutoff, the error-suppression benefits 
arising from the form of the filter function for UDD were reduced. 

4-2. Regimes of interest 

It is important to distinguish, at this point, between two salient regimes of operability 
in quantum systems. We note that a distinction between these regimes has not been 
formalized previously. 

(i) "High-fidelity" regime, occurring at times t << T 2 , where T 2 is the 1/e coherence 
time of the system. In this regime the accumulated error due to decoherence is 
small, but must be compared against, e.g. predicted fault-tolerance error thresholds 
of pth ~ 0.01% derived from quantum error correction. In quantum computing 
applications the maximum allowable error must not surpass pth, which is a much 
more stringent requirement than the typical experimental metric that r < T 2 . 

(ii) "Coherence time" regime, in which error probabilities of a few tens of percent 
are permissible, and T 2 is the only relevant metric. This may be appropriate 
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for, e.g. new quantum enabled sensing or imaging applications where large error 
probabilities provide contrast between different targets. 

4-2.1. The filter function in the high-fidelity regime The filter function interpretation 
of DD appropriate for the high-fidelity regime is shown graphically in Fig. [2] where the 
filter function of multiple n-pulse DD sequences is presented on a log-log plot for various 
values of n. In these plots we have normalized to a dimensionless angular frequency 



(to x r), as this is the only physical quantity that enters into Eq. 23. In all cases, the 
filter function is small at low frequencies and grows with frequency. Near ujpi it begins 
to oscillate rapidly. 

Examination of panels (a-c) of this figure show that the low-frequency rolloff in 
the stop-band varies dramatically between sequences, providing significant performance 
differences between these sequences for a given Sp(tS). Using the framework of filter- 
design we can analyze the filter functions for various sequences using new quantitative 
metrics. The low-frequency rolloff of PDD for these values of n is ~6 dB/octave while 
for the same n CPMG has a rolloff of ~18 dB/octave. This can be compared to the 
difference between the performance of single-pole and multi-pole filters in electrical 
engineering. As such, low-frequency noise is filtered much more efficiently using CPMG, 
where the filter function can be more than 100 dB lower than that for PDD given the 
same n. 

The value of up\ increases approximately linearly with n, indicating an increase in 
the spectral range efficiently filtered by the applied sequence. In neither case, however, 
does increasing pulse number change the low-frequency rolloff, i.e. the effective order 
of the filter. 

For both CPMG and PDD we observe an even-odd asymmetry (Fig. |2]i, e), but 
the pulse-number parity providing the best low-frequency rolloff differs between the 
sequences [521 [22] . In the case of PDD, even pulse numbers provide the worst suppression 
of low-frequency noise, comparable in form (~6 dB/octave) to that of FID. This is 
commensurate with the observation that given an even pulse number, in PDD we have 
an odd number of equal-length free-precession periods. As such, phase accumulation due 
to quasi-static noise is not fully compensated following the final free-precession period, 
making the system more susceptible to low-frequency noise. While this asymmetry has 
been recognized in the physics community, its origin is elucidated via the filter-design 
framework. 

The performance of the UDD sequence is quite distinct from the filtering 
characteristics of CPMG and PDD (Fig. [2|t). Most importantly, the UDD sequence 
provides a filter function whose rolloff increases in steepness with n. This is a 
manifestation of the analytic condition Uhrig originally imposed, in which the first n 
derivatives of F^ udd \ujt) = 0. As such, the addition of pulses in the UDD sequence is 
analogous to increasing the filter order in a standard filter design problem. By contrast, 
the low-frequency rolloff for CPMG and PDD is approximately constant with n, limiting 
their flexibility and performance. Additionally, even-odd parity differences in the UDD 
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filter functions are not visible in the slope of the low-frequency rolloff immediately below 

Differences in the ultimate error-suppressing capabilities of CPMG and UDD are 
clearly illustrated by a comparison of the filter functions for fixed n. In Fig. [3^, we 
plot the filter functions for n = 10 for CPMG and UDD, and the ratio of the filter 
functions in Fig. |3Jd. We clearly observe that for the majority of the spectral range of 
interest (below oupi), the filter function for UDD can be over 15 orders of magnitude 
smaller than that for CPMG, providing superior suppression of decoherence due to 
the associated components of Sp(u>). A broad spectral range of noise is "maximally" 
filtered by UDD before the onset of a sharp reduction of filter performance as u — > cJpi, 
an insight previously presented in various works. 

Based on these calculations, we see that for most instances, complex optimized 
sequence construction can provide significant benefits in noise suppression in the high- 
fidelity regime. UDD can outperform CPMG by many tens to hundreds of dB over 
experimentally relevant ranges of Sp(u), and given a sufficiently sharp high-frequency 
noise cutoff (corresponding to the sharp turn-on in the UDD filter function near Upi), 
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Figure 2. Graphic representation of the filter function for various pulse sequences as a 
function of dimensionless angular frequency (uxt). ( a_c ) Numerical filter function for 
various n and pulse sequences. Dashed grey line indicates the numerical value F — 1 
above which noise is not suppressed by the filter function, (d-f ) Filter functions for the 
same sequences and n — 2 — 6 demonstrating even-odd asymmetry in noise suppression 
of the various filter functions, as manifested in the slope of the low-frequency rolloff or 
the floor value in the case of UDD. 
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UDD generally gives near-optimal error suppressing performance in this regime. 

4-2.2. The filter function in the coherence time regime The tradeoff between spectral 
filtering range and filtering efficiency is also manifested in the observation that in most 
instances, for a given n, u^ 00 ^ 1 < ujp[ PMG ^ . We observe this in Fig. as a region in 
which CPMG outperforms UDD. Noise with significant power spectral density in the 
frequency range denoted by the dashed lines contributes to dephasing largely unimpeded 
as this region is close to uipi (Fig- 13^), or its effect may even be amplified by the presence 
of positive gain in the filter functions. As this noise contributes large dephasing-induced 
error probabilities, it predominantly impacts the measured value of T 2 in the system. 

Contributions to the system T 2 under dynamical decoupling sequence application 
may be clarified by wrapping the 1/u 2 term in Eq. [8] into the expression for the 
filter function [16J, yielding a modified filter function (ut) = F(cot)/co 2 . This 
construction emphasizes the parts of F(lot) near upi that provide the largest dephasing- 
error contributions and dominate T 2 . 

The modified filter functions for CPMG, PDD, and UDD are plotted on a linear 
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Figure 3. Comparison of the CPMG and UDD filter functions. (a)The filter functions 
calculated numerically for CPMG and UDD, with the value unity indicated by the 
dashed grey line, (b) The numerical value of the UDD filter function divided by that 
for CPMG. A value of one indicates equal filter functions, values greater than one 
correspond to CPMG outperforming UDD (red shading) , and vice versa for values less 
than one (green shading). Dashed lines in both panels indicate the frequency range 
near F = 1 where CPMG outperforms UDD, generally giving r £^ JDD ' 1 < r p^ CPMG% > 
in the presence of non-negligible noise spectral weight in this region. The benefits of 
UDD at the lowest frequencies are artificially limited as both filter functions saturate 
below ~ 10 -35 due to numerical limitations. 
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scale in Fig. [4], where pass-band ripple produces a peak in the modified filter function. 
This spectral peak serves as an approximate band-pass filter for noise dominating 
T2, with a minimum of ~20-26 dB of rejection outside the passband. Again, the 
pulse-sequences always serve as high-pass filters, but by focusing exclusively on the 
noise components that dominate the measured T 2 we find an approximate narrowband 
sampling of Sp(u). Note that for noise sufficiently dominated by low-frequency 
contributions, spectral components near ujpi will not necessarily dominate the measured 
T 2 for a given sequence. 

For CPMG and PDD we observe narrow spectral bands shifting towards higher 
frequencies with increasing n. This is derived from the fact that the use of equal 
interpulse free-precession times corresponds to selection of a fixed frequency component 
in the Fourier transform of the sampling function; in this case approximately a 
squarewave function. Low-frequency components are filtered out by the action of the 
pulses. This may be contrasted against the modified filter function for FID (Inset, 
Fig. Ek) that demonstrates significant contributions from noise at the lowest frequencies, 
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Figure 4. (a-c )Modified filter function, F^ 2 ^ (ujt) , for the pulse sequences and 
values of n studied in Fig. [2] on a linear scale. Dominant spectral contributions to the 
measured T2 appear as peaks in the modified filter functions. Inset a) Modified filter 
function for FID showing large weight for low-frequency noise on a semilog scale, with 
arbitrary units. Inset b) Demonstration of even-odd parity through the modified filter 
function of PDD on a semilog scale. 
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corresponding to extreme sensitivity to, e.g. shot-to-shot quasistatic fluctuations in 
an ensemble average measurement of error in a FID experiment. The low-frequency 
behavior of even-pulse-number-parity filter functions for PDD is demonstrated in the 
inset to panel (b). Like FID, low-frequency components are sampled broadband for even 
pulse numbers, equivalent to the discussion above about uncompensated quasistatic 
noise. Increasing the pulse number decreases the low-frequency constant value of 
the modified filter function. This is again expected as increasing n reduces the 
uncompensated phase accumulation time during the final free-precession period, thus 
reducing sensitivity to quasi-DC noise. These familiar behaviors are explained succinctly 
in the filter-design formalism. 

UDD again presents a departure from the observations for CPMG and PDD. In this 
case, increasing pulse number corresponds to a shift in spectral peak towards higher 
frequencies, but also a broadening of the spectral range contributing appreciably to 
T 2 . This arises through both a widening of the primary spectral peak with n, and 
also the emergence of higher-frequency oscillations above the primary peak. These 
features correspond to the slight increase in the standard filter function around ufi 
with n as depicted in Fig. ^p, and the fact that the sampling function of the UDD 
pulse sequence contains a variety of spectral components (corresponding to each of the 
different interpulse precession times). Broadband noise can therefore lead to a shortening 
of t( udd ^ relative to ^ CPMG ^ eV en if error rates at short times are lower for UDD than 
for CPMG. 

Interestingly, this form of the filter function has recently been employed in a pair of 
exciting experimental demonstrations. First, an experiment by Bylander and colleagues 
showed how the use of CPMG could permit tunable sampling of the noise power spectral 
density in a form of noise spectroscopy p2]. Second, Kotler et al demonstrated a 
"quantum lock-in amplifier" in which dynamical decoupling suppresses broadband noise 
but permits detection of a desired modulation signal chosen in the passband of a PDD 
sequence [61]. Our studies illustrate how, in the modified-filter function interpretation, 
periodic pulse sequences will provide the most narrow pass-band spectral response, 
while other sequences optimized for operation in the high-fidelity regime will sample 
noise broadband in contributing to T 2 . Thus a new "optimization" procedure will be 
needed for bandpass applications of DD sequences. 

5. Effects of Nonzero r w and other Timing Constraints 

Almost all theoretical studies of the effects of dynamical decoupling have been performed 
in the so-called "Bang-Bang" limit [21], in which control-pulse operations are assumed 
to be delta-functions, with pulse duration = 0. This assumption is convenient 
theoretically, but unphysical as any real system will have nonzero T n . The problem 
of accounting for realistic pulses [62] [63] has attracted considerable attention of late, 
and has yielded important insights. 

Returning to the filter function formalism we may make the lowest order 
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Figure 5. Modification of the filter function by inclusion of nonzero t w for UDD and 
CPMG with n = 7 and n — 8. In these panels different traces correspond to different 
values of the ratio of to the total free-precession time, r. 



approximation accounting for nonzero r n by modifying the sampling function y n (t) to 
incorporate a delay of length r n and value zero between values of ±1 corresponding 
to the free-precession periods [U [5] . This approximation assumes that environmentally 
induced phase accumulation is negligible during the application of a n pulse, consistent 
with many experimental observations in which Rabi decay times (the decay of driven 
oscillations) are much longer than FID times. When moving to the frequency domain, 
incorporating this delay results in a filter function for an arbitrary n-pulse sequence 



F(ut) 



+ 2£< 



-1)V 



cos (a;r 7r /2)|" 



(25) 



where Sjt' is the time of the center of the j th ttx pulse, and r' = r + nr n is the sum 
of the total free-precession time and 7r-pulse times. The underlying assumption of zero 
dephasing during applied pulses breaks down as r n — > r. In this limit, the system's 
evolution is dominated by periods of driven rotation rather than free-precession and 
evolution of f3{t) during extended pulses prevents the cancellation of phase accumulation 
in successive free-precession periods. 

The inclusion of such a delay in the sampling function, and the corresponding 
modification of the filter function is represented graphically in Fig. [5] for UDD and 
CPMG. Including a nonzero r n in this manner most significantly effects the UDD 
sequence in which the roll-off of the filter function is reduced by accounting for r n . 
As r n increases, the roll-off of the UDD filter function converges to approximately 
the same slope as that for CPMG, eliminating the benefits of the more complex 
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sequence construction, and identified for the first time here. For the regime studied 
in this figure, CPMG is minimally affected by the modification to the filter function. 
Here, tv is approximately six times smaller than the shortest interpulse free-precession 
period, t^^ mg \ but as t w becomes comparable to this quantity, ufi is shifted to lower 
frequencies (not shown). PDD behaves similarly to CPMG. While general modification 
of the filter function via inclusion of nonzero r n has been verified by experiment jU |5] , 
the new analyses of dynamical decoupling sequences as filters present significant new 
insight into the relative importance of r n for various sequences. 

A new analysis indicates that similar effects can be identified by reducing the 
precision with which pulse locations are specified (Fig. [6]). We may construct a 
filter function using instantaneous pulses but impose rounding on the fractional pulse 
locations used to produce the filter, effectively discretizing time and corresponding, 
e.g. to the effects of digital clocking. In particular, rounding of pulse timing accuracy 
to just one part in 10 7 can reduce the error suppression of the UDD filter function 
at a dimensionless angular frequency of unity by ~ 80 dB, with increasing losses in 
suppression at lower frequencies. This could correspond to the effects of a 10 GHz clock 
(lOOps clock period) on an experiment with duration 1 ms. Again, we observe that 
CPMG is nearly unaffected by the imposition of these constraints to the level of ~ 1%. 
The discretization of time discussed here may actually provide benefits in the analysis 
of filter design, borrowing from the control theory literature. 

Ultimately, this understanding of the filter function suggests that realistic hardware 
constraints may impose a nontrivial limit in the efficacy of optimized dynamical 
decoupling sequences such as UDD. Previous experiments have shown that such timing 
constraints could be ignored on the level of error rates ~ 10~ 3 , but may become 
appreciable as error rates are pushed near fault-tolerance limits [5]. 

6. Filter Optimization through Search Algorithms 

Each dynamical decoupling sequence discussed hereto had a specific filtering design 
objective; in the case of spin echo and CPMG that objective was to cancel slowly 
varying noise terms, while an n-pulse UDD sequence is constructed to cancel the first 
n orders in a Taylor expansion in time of any classical noise field These sequences 
have in common that they define a set of relative inter-pulse delays which is the same for 
arbitrary total sequence lengths. One might ask whether, by abandoning this philosophy, 
the filter design problem can be recast to find sequences that are either universally 
optimal or tailored to outperform all other sequences. We will see that it is possible to 
construct new suites of sequences, specifically tailored to a given noise environment in 
order to provide maximum error suppression. 
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Figure 6. The effect of imprecision in pulse location. The vector Sj is numerically 
rounded to the precision indicated in the legend. Maximum precision corresponds to 
double-precision integer representation. Here n — 10 and the main panel shows the 
filter function for UDD. Inset) Filter function for CPMG; effects of rounding are not 
visible up to 10 -2 



6.1. Locally Optimized Dynamical Decoupling (LODD) 

The first optimization approach involves designing a different noise filter for each value 
of the total sequence time, r, as this parameter sets the overall frequency characteristics 



of a given filter function (Eq. 23). By definition the best noise filter is that one which 
minimizes the coherence integral in EqJHj The corresponding analytical condition is 
that the derivative of the coherence integral with respect to each pulse location, Sj, 
must vanish, which leads to a set of coupled equations to solve. Since a solution to 
that set of equations is not possible for arbitrary Sp(u>) it is necessary to resort to 
numerical means. This philosophy was employed in jU |5] where the coherence integral 
was directly minimized for a given Sp(u) (as opposed to solving the set of coupled 
equations) via a multidimensional search algorithm. The resulting sequences were 
referred to as LODD sequences (Locally Optimized Dynamical Decoupling) in reference 
to the unique sequence obtained locally for each r. 

A set of such sequences, with n = 6, is shown in Fig. ([7]) where the open circles 
indicate the relative timings of the LODD sequence for each time indicated on the 
vertical axis. In this case the sequences were optimized for an ohmic spectrum S(u) = u 
with a sharp frequency cutoff at ujd- As with UDD, the LODD approach yields 
significant improvements when the noise spectrum is high-frequency dominated and 
has a sharp high-frequency cutoff jl]. Moreover LODD was shown to provide significant 
gains over UDD across a range of r under those circumstances. 

An advantage of the LODD approach is that it can be experimentally implemented 
using closed-loop feedback without any prior detailed knowledge of the power spectrum 
on the part of the operator; minimization of measured error is the only task to be 
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Figure 7. Relative pulse spacings of LODD (open circles) and OFDD (diamonds). 
The relative pulse spacing multiplied by the total sequence length, r, provides the 
pulse locations in absolute terms. 



carried out. This is accomplished by initiating a search for the optimized sequence with 
a random guess and measuring the resulting coherence. A new guess is then generated 
using a suitable search algorithm and the coherence remeasured. The process is repeated 
until the algorithm converges within the resolution afforded by the measurement fidelity. 



6.2. Optimized Noise Filtration Dynamical Decoupling (OFDD) 

The noise power spectrum enters directly into the numerical calculation of LODD 
sequences through the coherence integral, x{ T )- It is therefore reasonable to expect 
that the LODD sequences found would be very different for different noise spectra. One 
study has shown that locally optimized sequences can be found that are suitable for use 
under many circumstances [10] ■ Consider a new optimization objective of minimizing 



the area under the filter function, Eq. 23 up to some appropriate frequency cutoff, 



ujd- In an experimental system that cutoff will be set by a characteristic frequency of 
the system such as a spectral bandwidth of the noise, but need not be known a priori. 
Since the coherence integral is determined by the overlap between the filter function and 
noise spectrum, minimizing the area under the filter function will minimize the overlap 
without concern for detailed information about Sp(u). 

Again, a suite of pulse sequences is obtained, but this same set of sequences many 
be applied to a broad range of noise spectra. The pulse sequence at each total sequence 
length, r, is unique and predetermined numerically. This approach was referred to as 
optimized noise filtration dynamical decoupling (OFDD). In an experiment a feedback 
search algorithm need only perform a one-dimensional search through the predetermined 
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Figure 8. Comparison of filter functions for UDD (red), LODD (blue) and OFDD 
(black) sequences at r — 5 and t = 10. The dashed line corresponds to the imposed 
cutoff frequency, wjj. 

set of pulse sequences to find that one yielding the best coherence for a specific r, 
constituting a significantly reduced experimental-feedback burden. This then sets the 
optimum sequence for every other r. Figure [7] compares the pulse spacings for LODD 
optimized for the ohmic spectrum to that of the OFDD sequences in the same units. In 
the high-fidelity regime these sequences differ minimally. 

Reference [10] indicated that for certain noise environments the LODD and OFDD 
approaches can significantly improve on sequences in which the relative pulse spacings 
are the same for all sequence times, while doing no worse than other sequences in most 
environments. To gain further insight to why this is so it is instructive to compare 
the filter functions of locally optimized sequences directly. Figure [8] plots the filter 
functions for UDD (red), LODD (blue) and OFDD (black) sequences with n = 6 and 
at two different times r = 5 and r = 10 (units of 1/ojd). LODD was optimized for 
the same ohmic spectrum as in Fig. [7| Note that LODD and OFDD are qualitatively 
very similar. In both cases the filter functions are larger in magnitude than that due to 
UDD in the low-frequency regime, providing somewhat weaker noise suppression over 
this spectral range. Nonetheless these sequences perform better than UDD because they 
increase the slope of the low-frequency roll-off and hence suppress noise more strongly 
close to the high-frequency cutoff, i.e. that region where the filter function is usually 
largest, and where for an Ohmic spectrum the noise is strongest. 
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6.3. Bandwidth Adapted Dynamical Decoupling (BADD) 

In experimental settings control technologies impose constraints on DD often neglected 
in theoretical studies, such as imperfect detection, pulse errors, finite energy constraints, 
and finite switching time of the control modulation. Incorporating such control 
constraints into an optimization routine can produce novel sequences which improve 
dephasing supression. For instance, the LODD sequence construction described above 
implicitly accounts for r n by simply finding the suite of pulse sequences minimizing 
x(t'). A more explicit accounting for timing constraint is addressed by Khodjasteh et 
al., in the formulation of Bandwidth-Adapted Dynamical Decoupling (BADD) |41j . 

There, Khodjasteh et al. [JT] re-evaluated the numerical DD optimization problem 
in light of finite timing resources; they imposed the constraint that the shortest inter- 
pulse delay, r min , in any DD sequence must obey: 

T~min ^ T switch-, (^6) 

where T SW i tc h is a minimum switching time allowed by the control technology. They 
then performed numerical searches for optimum pulse sequences in which both r m i n 
and the total sequence time r are constrained, but n is allowed to vary. The resulting 
sequences were referred to as bandwidth- adapted DD (BADD) sequences. With these 
timing constraints UDD sequences are much more severely limited in the number of 
allowed pulses for a given r than the BADD approach. 

The authors compare the decoupling error for UDD and BADD sequences in a spin- 
boson dephasing model with a supra-Ohmic power spectrum Sp{u) = aa/exp (—uj/uj c ), 
where a ~ 1.14 x 1CT 26 s 2 and u c ~ 3 rad/ps. There a timing constraint is of T min ^0.1 
ps is imposed by the need to avoid unwanted excitation of higher energy levels. Then for 
a total sequence time of r ^ 10 2 ps BADD and LODD sequences allow up to n = 100 
while UDD can only accomodate n = 20 commensurate with the constraint. The authors 
observe that the minimum achievable error is only very weakly dependent on the total 
sequence time for the BADD sequences. This is an indication that when operating in a 
regime of 'fast-control' DD performance should be measured in terms of T SW i tc h rather 
than r, to allow fair comparisons. UDD on the other hand performs significantly poorer 
due to the n = 20 limit. 

7. Conclusions 

In this manuscript we have discussed the construction of dynamical decoupling pulse 
sequences for the efficient suppression of decoherence as a problem in filter design, 
analogous to tasks in digital signal processing or analog electronics. This approach 
provides a practical perspective on the functionality of various dynamical decoupling 
pulse sequences, and helps to elucidate the relative performance advantages and 
limitations of certain optimized constructions. This work has provided a comprehensive 
review of recent results in the field, but has also yielded new results derived explicitly 
from the use of a filter-design analytical framework. 
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We have analyzed the performance of the leading pulse sequence formulations 
including PDD, CPMG, and UDD, and reviewed new optimization techniques such 
as LODD, OFDD, and BADD entirely in the context of filter design. A quantitative 
analysis derived from classical control theory forms the theoretical foundation for this 
treatment, and is based upon the reduction of dephasing-suppression to a problem 
in linear control theory. Our analyses have yielded new insights into how these 
sequences provide robust performance in particular noise environments, and they allow 
an experimentalist to accurately predict error suppression performance and select an 
optimized sequence given ambient conditions. Specifically, we have identified the filter 
roll-off and ripple as limiting performance for various DD sequences, and studied how 
changes in pulse sequence parameters such as the parity of n impact the roll-off. Our 
discussion also included the consideration of realistic experimental constraints such as 
nonzero pulse durations and novel conclusions about timing constraints, describing 
potential impacts on filter performance. We note that similar general insights - 
without the specific use of the filter-design framework - have recently appeared in the 
literature [39] 

We expect continued exciting developments in dynamical decoupling sequence 
construction as the community begins to leverage the decades of insight provided by the 
field of signal processing for robust filter design. In principle, insights coming directly 
from the classical control community may be leveraged for the creation of new sequences 
in the future, but producing useful sequences in this way remains the subject of future 
work. 

In the perspective presented in this paper we ignored concatenated dynamical 
decoupling as discussed in (3U1 EU ESI El ES] • One might be tempted to assume that using 
local optimization techniques one could always find pulse sequences that can outperform 
any CDD scheme for a given number of pulses. However, it is important to remember 
that with increased pulse numbers the optimization burden also rapidly increases. 
Therefore, new efforts focused on the creation of optimized pulse sequences for the 
robust preservation of arbitrary qubit states will likely incorporate optimized sequence 
construction for each level of a concatenated sequence as suggested in Refs. [42, 43]. 
The combination of these two approaches promises to provide vital capabilities for the 
implementation of error-resistent quantum memory. 

Future work will also include the transfer of insights gained from dynamical 
decoupling sequence construction for implementation of the identity operator towards 
dynamical error suppression during the application of nontrivial logic - dynamically 
corrected gates [HZl ESI ESI ED] • Combined with new enhanced experimental capabilities, 
we believe efficient dynamical error suppression will become a foundational underlying 
functionality in many quantum coherent technologies. 
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